Pharmacophore Variants of the Macrocyclic Peptide Triazole Inactivator of HIV-1 Env.

Previously we established a family of macrocyclic peptide triazoles (cPTs) that inactivate the Env protein complex of HIV-1, and identified the pharmacophore that engages Env's receptor binding pocket. Here, we examined the hypothesis that the side chains of both components of the triazole Pro - Trp segment of cPT pharmacophore work in tandem to make intimate contacts with two proximal subsites of the overall CD4 binding site of gp120 to stabilize binding and function. Variations of the triazole Pro R group, which previously had been significantly optimized, led to identification of a variant MG-II-20 that contains a pyrazole substitution. MG-II-20 has improved functional properties over previously examined variants, with Kd for gp120 in the nM range. In contrast, new variants of the Trp indole side chain, with either methyl- or bromo- components appended, had disruptive effects on gp120 binding, reflecting the sensitivity of function to changes in this component of the encounter complex. Plausible in silico models of cPT:gp120 complex structures were obtained that are consistent with the overall hypothesisof occupancy by the triazole Pro and Trp side chains, respectively, into the ß20/21 and Phe43 sub-cavities. The overall results strengthen the definition of the cPT-Env inactivator binding site and provide a new lead composition (MG-II-20) as well as structure-function findings to guide future HIV-1 Env inactivator design.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels.

The recent determination of cryo-EM structures of voltage-gated sodium (Nav) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Nav pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Nav channels based on homology modeling of the cryo-EM structure of the human Nav1.4 channel and, in addition, on the recently resolved configuration for Nav1.2. In particular, single Na+ permeation events during standard MD runs suggest that the ion resides in the inner part of the Nav selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na+ translocation through the SF of the homology-based Nav1.2 model and the cryo-EM Nav1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na+ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Nav cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Nav channels.

Structural and Functional Characterization of Indane-Core CD4-Mimetic Compounds Substituted with Heterocyclic Amines.

The human immunodeficiency virus (HIV-1) envelope glycoprotein (Env) trimer on the virion surface interacts with the host receptors, CD4 and CCR5/CXCR4, to mediate virus entry into the target cell. CD4-mimetic compounds (CD4mcs) bind the gp120 Env, block CD4 binding, and inactivate Env. Previous studies suggested that a C(5)-methylamino methyl moiety on a lead CD4mc, BNM-III-170, contributed to its antiviral potency. By replacing the C(5) chain with differentially substituted pyrrolidine, piperidine, and piperazine ring systems, guided by structural and computational analyses, we found that the 5-position of BNM-III-170 is remarkably tolerant of a variety of ring sizes and substitutions, both in regard to antiviral activity and sensitization to humoral responses. Crystallographic analyses of representative analogues from the pyrrolidine series revealed the potential for 5-substituents to hydrogen bond with gp120 Env residue Thr 283. Further optimization of these interactions holds promise for the development of CD4mcs with greater potency.

Entropic Overcompensation of the N501Y Mutation on SARS-CoV-2 S Binding to ACE2.

Recent experimental work has shown that the N501Y mutation in the SARS-CoV-2 S glycoprotein's receptor binding domain (RBD) increases binding affinity to the angiotensin-converting enzyme 2 (ACE2), primarily by overcompensating for a less favorable enthalpy of binding by greatly reducing the entropic penalty for complex formation, but the basis for this entropic overcompensation is not clear [Prévost et al. J. Biol. Chem.2021, 297, 101151]. We use all-atom molecular dynamics simulations and free-energy calculations to qualitatively assess the impact of the N501Y mutation on the enthalpy and entropy of binding of RBD to ACE2. Our calculations correctly predict that N501Y causes a less favorable enthalpy of binding to ACE2 relative to the original strain. Furthermore, we show that this is overcompensated for by a more entropically favorable increase in large-scale quaternary flexibility and intraprotein root mean square fluctuations of residue positions upon binding in both RBD and ACE2. The enhanced quaternary flexibility stems from N501Y's ability to remodel the inter-residue interactions between the two proteins away from interactions central to the epitope and toward more peripheral interactions. These findings suggest that an important factor in determining protein-protein binding affinity is the degree to which fluctuations are distributed throughout the complex and that residue mutations that may seem to result in weaker interactions than their wild-type counterparts may yet result in increased binding affinity thanks to their ability to suppress unfavorable entropy changes upon binding.

Entropic overcompensation of the N501Y mutation on SARS-CoV-2 S binding to ACE2.

Recent experimental work has shown that the N501Y mutation in the SARS-CoV-2 S glycoprotein's receptor binding domain (RBD) increases binding affinity to the angiotensin-converting enzyme 2 (ACE2), primarily by overcompensating for a less favorable enthalpy of binding by a greatly reducing the entropic penalty for complex formation, but the basis for this entropic overcompensation is not clear [Prévost et al., J. Biol. Chem . (2021) 297;101151]. We use all-atom molecular dynamics simulations and free-energy calculations to qualitatively assess the impact of the N501Y mutation on enthalpy and entropy of binding of RBD to ACE2. Our calculations correctly predict that N501Y causes a less favorable enthalpy of binding to ACE2 relative to the original strain. Further, we show that this is overcompensated for by a more entropically favorable increase in large-scale quaternary flexibility and intra-protein root-mean squared fluctuations of residue positions upon binding in both RBD and ACE2. The enhanced quaternary flexibility stems from N501Y's ability to remodel the interresidue interactions between the two proteins away from interactions central to the epitope and toward more peripheral interactions. These findings suggest that an important factor in determining protein-protein binding affinity is the degree to which fluctuations are distributed throughout the complex, and that residue mutations that may seem to result in weaker interactions than their wild-type counterparts may yet result increased binding affinity thanks to their ability to suppress unfavorable entropy changes upon binding.

Target identification for repurposed drugs active against SARS-CoV-2 via high-throughput inverse docking.

Screening already approved drugs for activity against a novel pathogen can be an important part of global rapid-response strategies in pandemics. Such high-throughput repurposing screens have already identified several existing drugs with potential to combat SARS-CoV-2. However, moving these hits forward for possible development into drugs specifically against this pathogen requires unambiguous identification of their corresponding targets, something the high-throughput screens are not typically designed to reveal. We present here a new computational inverse-docking protocol that uses all-atom protein structures and a combination of docking methods to rank-order targets for each of several existing drugs for which a plurality of recent high-throughput screens detected anti-SARS-CoV-2 activity. We demonstrate validation of this method with known drug-target pairs, including both non-antiviral and antiviral compounds. We subjected 152 distinct drugs potentially suitable for repurposing to the inverse docking procedure. The most common preferential targets were the human enzymes TMPRSS2 and PIKfyve, followed by the viral enzymes Helicase and PLpro. All compounds that selected TMPRSS2 are known serine protease inhibitors, and those that selected PIKfyve are known tyrosine kinase inhibitors. Detailed structural analysis of the docking poses revealed important insights into why these selections arose, and could potentially lead to more rational design of new drugs against these targets.

Identification of gp120 Residue His105 as a Novel Target for HIV-1 Neutralization by Small-Molecule CD4-Mimics.

The design and synthesis of butyl chain derivatives at the indane ring 3-position of our lead CD4-mimetic compound BNM-III-170 that inhibits human immunodeficiency virus (HIV-1) infection are reported. Optimization efforts were guided by crystallographic and computational analysis of the small-molecule ligands of the Phe43 cavity of the envelope glycoprotein gp120. Biological evaluation of 11-21 revealed that members of this series of CD4-mimetic compounds are able to inhibit HIV-1 viral entry into target cells more potently and with greater breadth compared to BNM-III-170. Crystallographic analysis of the binding pocket of 14, 16, and 17 revealed a novel hydrogen bonding interaction between His105 and a primary hydroxyl group on the butyl side chain. Further optimization of this interaction with the His105 residue holds the promise of more potent CD4-mimetic compounds.

Atomistic simulation of volumetric properties of epoxy networks: effect of monomer length.

Properties of epoxy thermosets can be varied broadly to suit design requirements by altering the chemistry of the component agents. Atomistically-detailed molecular dynamics simulations are well-suited for molecular insight into the structure-property relationship for a rational tailoring of the chemistry. Since the macroscopic properties of interest for applications emerge hierarchically from molecular-scale chemical interactions, seamless integration of experiment, computation, and theory is of great interest. Recently, a Specific Volume-Cooling Rate analysis protocol was successfully developed to quantitatively compare the volumetric properties of an epoxy network model with experimental results in the literature, in spite of the nine orders of magnitude mismatch in the accessible time-scales. Here, we extend the application of the method for two epoxy networks in the same class of chemistry but whose monomers have a higher number of repeating units compared to the previous one for validating the generality of our approach. We observed that atomistic simulations are able to predict the experimental temperature trend of the specific volume within 0.4% for both these networks. Using the William-Landel-Ferry equation to account for rate differences, we also see good agreement between the computational and experimental values of the glass transition temperature.

Impact of temperature on the affinity of SARS-CoV-2 Spike glycoprotein for host ACE2.

The seasonal nature of outbreaks of respiratory viral infections with increased transmission during low temperatures has been well established. Accordingly, temperature has been suggested to play a role on the viability and transmissibility of SARS-CoV-2, the virus responsible for the COVID-19 pandemic. The receptor-binding domain (RBD) of the Spike glycoprotein is known to bind to its host receptor angiotensin-converting enzyme 2 (ACE2) to initiate viral fusion. Using biochemical, biophysical, and functional assays to dissect the effect of temperature on the receptor-Spike interaction, we observed a significant and stepwise increase in RBD-ACE2 affinity at low temperatures, resulting in slower dissociation kinetics. This translated into enhanced interaction of the full Spike glycoprotein with the ACE2 receptor and higher viral attachment at low temperatures. Interestingly, the RBD N501Y mutation, present in emerging variants of concern (VOCs) that are fueling the pandemic worldwide (including the B.1.1.7 (α) lineage), bypassed this requirement. This data suggests that the acquisition of N501Y reflects an adaptation to warmer climates, a hypothesis that remains to be tested.

Roles of the Coupling Agent and Surfactant in Droplet Structure in Sizing Emulsions: A Molecular Dynamics Simulations Study.

Sizing emulsions used as glass fiber surface treatments in composites manufacturing are aqueous suspensions of hydrophobic film formers, surface coupling agents, and surfactants. We employ all-atom molecular dynamics simulations to characterize droplet structures in several aqueous blends of the film-former diglycidyl ether of bisphenol A, coupling agent glycidoxypropyl trimethoxysilane, and a triblock copolymer surfactant (Pluronic L35 PEO/PPO copolymer). We show that the quasi-equilibrium states of emulsion droplets are invariant to different initial configurations. We examine the role of the surfactant in determining coupling agent partitioning between the droplet shell and corona and coupling agent cluster size distributions. This work takes a step toward systematic understanding of the sizing chemistry to optimize the interface between the glass and the resin in commercially relevant composites.

Peptide Triazole Thiol Irreversibly Inactivates Metastable HIV-1 Env by Accessing Conformational Triggers Intrinsic to Virus-Cell Entry.

KR13, a peptide triazole thiol previously established to inhibit HIV-1 infection and cause virus lysis, was evaluated by flow cytometry against JRFL Env-presenting cells to characterize induced Env and membrane transformations leading to irreversible inactivation. Transiently transfected HEK293T cells were preloaded with calcein dye, treated with KR13 or its thiol-blocked analogue KR13b, fixed, and stained for gp120 (35O22), MPER (10E8), 6-helix-bundle (NC-1), immunodominant loop (50-69), and fusion peptide (VRC34.01). KR13 induced dose-dependent transformations of Env and membrane characterized by transient poration, MPER exposure, and 6-helix-bundle formation (analogous to native fusion events), but also reduced immunodominant loop and fusion peptide exposure. Using a fusion peptide mutant (V504E), we found that KR13 transformation does not require functional fusion peptide for poration. In contrast, simultaneous treatment with fusion inhibitor T20 alongside KR13 prevented membrane poration and MPER exposure, showing that these events require 6-helix-bundle formation. Based on these results, we formulated a model for PTT-induced Env transformation portraying how, in the absence of CD4/co-receptor signaling, PTT may provide alternate means of perturbing the metastable Env-membrane complex, and inducing fusion-like transformation. In turn, the results show that such transformations are intrinsic to Env and can be diverted for irreversible inactivation of the protein complex.

Impact of temperature on the affinity of SARS-CoV-2 Spike for ACE2.

The seasonal nature in the outbreaks of respiratory viral infections with increased transmission during low temperatures has been well established. The current COVID-19 pandemic makes no exception, and temperature has been suggested to play a role on the viability and transmissibility of SARS-CoV-2. The receptor binding domain (RBD) of the Spike glycoprotein binds to the angiotensin-converting enzyme 2 (ACE2) to initiate viral fusion. Studying the effect of temperature on the receptor-Spike interaction, we observed a significant and stepwise increase in RBD-ACE2 affinity at low temperatures, resulting in slower dissociation kinetics. This translated into enhanced interaction of the full Spike to ACE2 receptor and higher viral attachment at low temperatures. Interestingly, the RBD N501Y mutation, present in emerging variants of concern (VOCs) that are fueling the pandemic worldwide, bypassed this requirement. This data suggests that the acquisition of N501Y reflects an adaptation to warmer climates, a hypothesis that remains to be tested.

Optimizing String Method's Reproducibility Using Generalized Solute Tempering Replica Exchange.

Obtaining accurate and reproducible free energies from molecular simulations is somewhat tricky due to incomplete knowledge of crucial slow degrees of freedom leading to hidden barriers that can stymie sampling. Employing a sufficiently large number of collective variables (CV) and ensuring ergodic sampling in orthogonal CV space, perhaps via tempering methods, can reduce these issues to some extent. For complex systems with high-dimensional free energy landscapes, both these approaches become computationally expensive. For high-dimensional landscapes, efficient exploration can be enabled by using temperature-accelerated MD (TAMD) and identification and characterization of minimum free energy pathways connecting minima can be found by using the string method (SM). Both TAMD and SM use mean-force estimates from finite MD simulations and are thus susceptible to sampling restrictions from hidden variables. A recent development in parallel tempering methods, "generalized replica exchange solute tempering" (gREST), can enhance sampling at a reasonable computational cost with its flexibility to target very specific "solutes" which can include arbitrary independent variables. Considering the advantages of both methods, we implement gREST-enabled TAMD and SM. By considering two different collective variable representations of the pentapeptide neurotransmitter met-enkephalin, we show that both gREST-enabled TAMD and SM yield more accurate and reproducible free energy predictions than TAMD and SM alone. Given the moderate computational cost of gREST compared with other replica-exchange methods, gREST-enabled SM represents a more attractive method for characterizing free energy minima and pathways among them for a large variety of systems.

HIV-1 Env-Dependent Cell Killing by Bifunctional Small-Molecule/Peptide Conjugates.

A strategy has been established for the synthesis of a family of bifunctional HIV-1 inhibitor covalent conjugates with the potential to bind simultaneously to both the gp120 and gp41 subunits of the HIV-1 envelope glycoprotein trimeric complex (Env). One component of the conjugates is derived from BNM-III-170, a small-molecule CD4 mimic that binds to gp120. The second component, comprised of the peptide DKWASLWNW ("Trp3"), was derived from the N-terminus of the HIV-1 gp41 Membrane Proximal External Region (MPER) and found previously to bind to the gp41 subunit of Env. The resulting bifunctional conjugates were shown to inhibit virus cell infection with low micromolar potency and to induce lysis of the HIV-1 virion. Crucially, virolysis was found to be dependent on the covalent linkage of the BNM-III-170 and Trp3 domains, as coadministration of a mixture of the un-cross-linked components proved to be nonlytic. However, a significant magnitude of lytic activity was observed in Env-negative and other control pseudoviruses, suggesting parallel mechanisms of action of the conjugates involving Env interaction and direct membrane disruption. Computational modeling suggested strong membrane-binding activity of BNM-III-170, which may underly the nonspecific virolytic effects of the conjugates. To investigate the scope of the membrane effect, cell-based cytotoxicity and membrane permeability assays were performed employing flow cytometry. Here, we observed a dose-dependent and specific cytotoxic effect on HIV-1 Env-expressing cells by the small-molecule bifunctional inhibitor. Most importantly, Env-negative cells were not susceptible to the cytotoxic effect upon exposure to this construct at concentrations where cell-killing effects were observed for Env-positive cells. Computational structural modeling supports a mechanism in which the bifunctional inhibitors bind to the gp120 and gp41 subunits in tandem in open-state Env trimers and induce relative motion of the gp120 subunits consistent with models of Env inactivation. This observation supports the idea that the cell-killing effect of the small-molecule bifunctional inhibitor is due to specific Env conformational triggering. This work lays important groundwork to advance a small-molecule bifunctional inhibitor approach for eliminating Env-expressing infected cells and the eradication of HIV-1.

Interaction of Human ACE2 to Membrane-Bound SARS-CoV-1 and SARS-CoV-2 S Glycoproteins.

Severe acute respiratory syndrome virus 2 (SARS-CoV-2) is responsible for the current global coronavirus disease 2019 (COVID-19) pandemic, infecting millions of people and causing hundreds of thousands of deaths. The viral entry of SARS-CoV-2 depends on an interaction between the receptor-binding domain of its trimeric spike glycoprotein and the human angiotensin-converting enzyme 2 (ACE2) receptor. A better understanding of the spike/ACE2 interaction is still required to design anti-SARS-CoV-2 therapeutics. Here, we investigated the degree of cooperativity of ACE2 within both the SARS-CoV-2 and the closely related SARS-CoV-1 membrane-bound S glycoproteins. We show that there exist differential inter-protomer conformational transitions between both spike trimers. Interestingly, the SARS-CoV-2 spike exhibits a positive cooperativity for monomeric soluble ACE2 binding when compared to the SARS-CoV-1 spike, which might have more structural restraints. Our findings can be of importance in the development of therapeutics that block the spike/ACE2 interaction.

Identification of a glycan cluster in gp120 essential for irreversible HIV-1 lytic inactivation by a lectin-based recombinantly engineered protein conjugate.

We previously discovered a class of recombinant lectin conjugates, denoted lectin DLIs ('dual-acting lytic inhibitors') that bind to the HIV-1 envelope (Env) protein trimer and cause both lytic inactivation of HIV-1 virions and cytotoxicity of Env-expressing cells. To facilitate mechanistic investigation of DLI function, we derived the simplified prototype microvirin (MVN)-DLI, containing an MVN domain that binds high-mannose glycans in Env, connected to a DKWASLWNW sequence (denoted 'Trp3') derived from the membrane-associated region of gp41. The relatively much stronger affinity of the lectin component than Trp3 argues that the lectin functions to capture Env to enable Trp3 engagement and consequent Env membrane disruption and virolysis. The relatively simplified engagement pattern of MVN with Env opened up the opportunity, pursued here, to use recombinant glycan knockout gp120 variants to identify the precise Env binding site for MVN that drives DLI engagement and lysis. Using mutagenesis combined with a series of biophysical and virological experiments, we identified a restricted set of residues, N262, N332 and N448, all localized in a cluster on the outer domain of gp120, as the essential epitope for MVN binding. By generating these mutations in the corresponding HIV-1 virus, we established that the engagement of this glycan cluster with the lectin domain of MVN*-DLI is the trigger for DLI-derived virus and cell inactivation. Beyond defining the initial encounter step for lytic inactivation, this study provides a guide to further elucidate DLI mechanism, including the stoichiometry of Env trimer required for function, and downstream DLI optimization.

Roles of variable linker length in dual acting virucidal entry inhibitors on HIV-1 potency via on-the-fly free energy molecular simulations.

The Dual-Acting Virolytic Entry Inhibitors, or DAVEI's, are a class of recombinant chimera fusion proteins consisting of a lectin, a flexible polypeptide linker, and a fragment of the membrane-proximal external region (MPER) of HIV-1 gp41. DAVEIs trigger virolysis of HIV-1 virions through interactions with the trimeric envelope glycoprotein complex (Env), though the details of these interactions are not fully determined as yet. The purpose of this work was to use structural modeling to rationalize a dependence of DAVEI potency on the molecular length of the linker connecting the two components. We used temperature accelerated molecular dynamics and on-the-fly parameterization to compute free energy versus end-to-end distance for two different linker lengths, DAVEI L0 (His6 ) and DAVEI L2 ([Gly4 Ser]2 His6 ). Additionally, an envelope model was created based on a cryo-electron microscopy-derived structure of a cleaved, soluble Env construct, with high-mannose glycans added which served as putative docking locations for the lectin, along with MPER added that served as a putative docking location for the MPER region of DAVEI (MPERDAVEI ). Using MD simulation, distances between the lectin C-terminus and Env gp41 MPER were measured. We determined that none of the glycans were close enough to gp41 MPER to allow DAVEI L0 to function, while one, N448, will allow DAVEI L2 to function. These findings are consistent with the previously determined dependence of lytic function on DAVEI linker lengths. This supports the hypothesis that DAVEI's engage Env at both glycans and the Env MPER in causing membrane poration and lysis.

The HIV-1 Env gp120 Inner Domain Shapes the Phe43 Cavity and the CD4 Binding Site.

The HIV-1 envelope glycoproteins (Env) undergo conformational changes upon interaction of the gp120 exterior glycoprotein with the CD4 receptor. The gp120 inner domain topological layers facilitate the transition of Env to the CD4-bound conformation. CD4 engages gp120 by introducing its phenylalanine 43 (Phe43) in a cavity ("the Phe43 cavity") located at the interface between the inner and outer gp120 domains. Small CD4-mimetic compounds (CD4mc) can bind within the Phe43 cavity and trigger conformational changes similar to those induced by CD4. Crystal structures of CD4mc in complex with a modified CRF01_AE gp120 core revealed the importance of these gp120 inner domain layers in stabilizing the Phe43 cavity and shaping the CD4 binding site. Our studies reveal a complex interplay between the gp120 inner domain and the Phe43 cavity and generate useful information for the development of more-potent CD4mc.IMPORTANCE The Phe43 cavity of HIV-1 envelope glycoproteins (Env) is an attractive druggable target. New promising compounds, including small CD4 mimetics (CD4mc), were shown to insert deeply into this cavity. Here, we identify a new network of residues that helps to shape this highly conserved CD4 binding pocket and characterize the structural determinants responsible for Env sensitivity to small CD4 mimetics.

Optimization of Small Molecules That Sensitize HIV-1 Infected Cells to Antibody-Dependent Cellular Cytotoxicity.

With approximately 37 million people living with HIV worldwide and an estimated 2 million new infections reported each year, the need to derive novel strategies aimed at eradicating HIV-1 infection remains a critical worldwide challenge. One potential strategy would involve eliminating infected cells via antibody-dependent cellular cytotoxicity (ADCC). HIV-1 has evolved sophisticated mechanisms to conceal epitopes located in its envelope glycoprotein (Env) that are recognized by ADCC-mediating antibodies present in sera from HIV-1 infected individuals. Our aim is to circumvent this evasion via the development of small molecules that expose relevant anti-Env epitopes and sensitize HIV-1 infected cells to ADCC. Rapid elaboration of an initial screening hit using parallel synthesis and structure-based optimization has led to the development of potent small molecules that elicit this humoral response. Efforts to increase the ADCC activity of this class of small molecules with the aim of increasing their therapeutic potential was based on our recent cocrystal structures with gp120 core.